Inelastic cotunneling through a long diffusive wire 
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We show that electron transport through a long multichannel wire, connected to leads by tunnel 
junctions, at low temperatures T and voltages V is dominated by inelastic cotunnelling. This 
mechanism results in experimentally observed power-law dependence of conductance on T and V, 
in the diffusive regime where usual Coulomb anomaly theory leads to exponentially low conductance. 
The power-law exponent a* is proportional to the distance between contacts L. 

PACS numbers: 73.63.-b, 73.23.Hk 



Electronic transport through nanowires was intensively 
studied by many groups in the past years. In particular, 
conductance of mult i- channel diffusive nanowires with 
relatively poor contacts to metal terminals was measured, 
cf. e.g. The Coulomb phenomena play an important 
role in transport, provided the contacts between the wire 
and the leads are weak. The mechanism of the Coulomb 
blockade, as well as the Coulomb anomaly due to tunnel- 
ing spreading of charge, are presently well understood. 
In the ballistic regime (at relatively high temperature T 
and/or bias voltage V) the Coulomb effects lead to the 
power-law temperature and voltage dependence of the 
conductance: 

G = dI/dV (xV"' (lowT), GcxT" (high T), (1) 

characteristic for Luttinger Liquid, while at low T and 
F - in diffusive regime - an exponential dependence (see 
(|4]) below) should be observed. The puzzle is that the 
power law ([1]) is found in almost all experiments, even in 
those where the conditions for the diffusive regime seem 
to be fulfilled. 

The existing theories (see, e.g., 01 Hi 01) considered 
the Coulomb effects at each of contacts separately. How- 
ever, if both contacts are taken into account simultane- 
ously, then some analog of cotunneling becomes possible 
and at low temperatures this mechanism should dom- 
inate. The standard theory of cotunneling deals with 
small grains or quantum dots, while a long wire is an ex- 
tended object: internal dynamics of charge within it may 
be important. In the present letter we develop a theory 
for such an extended cotunneling and show that in the 
diffusive regime the resulting cotunneling conductance 
still obeys the law ([T]), though with different exponent 
a* , depending on the separation L between the contacts. 

Consider a multichannel metallic wire (it may be a 
multiwall nanotube) of length Lq and diameter a. The 
wire is connected to massive metallic leads through two 
weak tunnel contacts A and B with identical dimension- 
less conductances 5^1 (see FiglT]). A voltage, apphed 
between the contacts is V. The classic dimensionless re- 
sistance of the piece of wire between the leads is assumed 
to be not very small: R{L)/{h/e'^) = L/^ > 1, where 



u 

a d 

1 



FIG. 1: Electrons tunnel between a wire of length Lq and 
diameter a and two leads A and B, placed symmetrically with 
respect to the center of the wire, at distance L from each 
other. 



C ^ Nchl is the localization length, I is the mean free 
path, and iVch 3> 1 is the number of channels. 

In this paper the relevant energy scales will be assumed 
so low, that the motion of electrons in the wire is diffu- 
sive. On the other hand, we will neglect the localization 
effects. As long as usual conductivity in a wire is con- 
cerned, the condition of " no localization" reads 



vf/K 



(2) 



It is not evident that inequlity ^ is in fact necessary 
when the under-the-barrier spreading process is consid- 
ered; however, it is certainly the sufficient one, and we 
will assume it is fulfilled below. This requirement is con- 
sistent at A'ch ^ 1 with the diffusive dynamics of charge 
spreading. 

If the temperature T is not very low, the diffusive 
transport between the two leads proceeds in a "single- 
particle mode": At first one electron (one hole) tunnels 
into the wire from one of the leads and is accommodated 
in the wire, then one hole (one electron) tunnels from an- 
other lead. Because of the (thermoactivated) tunneling 
character of the accommodation process, the correspond- 
ing conductance is exponentially suppressed 



G 



(1) 

AB 



gcxp{^Si{T,V)}. 



(3) 



At temperatures T :$> Tc = EcH) accommodation 
proceeds according to the semiclassic scenario [21 (the 
Coulomb zero-bias-anomaly regime, see also [3,11, 0]), 
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and the accommodation action 



SiiT,V) 



EciO/ieV), for » /B^(OT- 



(4) 



Here Ec(x) = \n.{x / a) / ex is the charging energy of a 
piece of wire of length x 3> a. 

At T < ri^^ the single-particle accomodation pro- 
ceeds according the "orthodox" Coulomb blockade sce- 
nario (see, e.g., 0): 



Si{T, V)^Ec{Lo)/T. 



(5) 



The abovementioned independent single particle pro- 
cesses should be less effective than some correlated cotun- 
neling process, in which any charged states of the wire 
would only enter as virtual intermediate states. The the- 
ory of such processes is well developed for transitions via 
small grains, where the intergrain charge transfer pro- 
cesses are the bottlenecks for the trarisport, while the 
intragrain charge transfer is easy (see [8|). In our case, 
however, the charge spreading within the wire is a crucial 
ingredient of the process, so that the standard perturba- 
tional description of the cotunneling is inapplicable. 

In the present Letter we propose a modification of 
the approach which allows for description of charge 
spreading effects under the two-particle cotunneling con- 
ditions. Our main result reads as follows: 
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(2) 
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max{LT/C, {eV)} 
Ec{L) 



.(6) 



In the case L ^ Lq the crossover from one-particle tunnel- 
ing (in the Coulomb blockade mode) to the two-particle 
one takes place at 



y(2) 



ln(io/a) 



^^1 In 



ln(l/.g) 



< 



(7) 



The last inequality becomes strong for very low conduc- 
tance of contacts, 3 <C 1; in this situation a sequence of 
crossovers may be seen with the temperature decrease: 
from the Coulomb anomaly mode (|4]) to the Coulomb 
blockade mode (O at T = Tc^\ and then to the inelastic 
cotunneling regime (O at T = Tc^\ 

In the case L Lq the Coulomb blockade regime is 
absent, and the crossover from the Coulomb anomaly to 
the inelastic cotunneling takes place at 



y(2) 



ln(^/a) 



£L2 



In 



+ 4ln(l/ff) 



(8) 



In the nonlinear regime the crossover between the single- 
particle Coulomb anomaly and the inelastic cotunneling 
takes place at 



eV ~ (e1/)c 



eL 



ln(C/a). 



(9) 



Thus, at low enough temperatures, the cotunneling sce- 
nario always dominates. On the other hand, the condi- 
tion ^ of "no localization" should also be fulfilled for 
applicability of the formula ^ . The necessary tempera- 
ture range only exists if 




ln(C/a) 



N, 



ch 



1/2 



(10) 



i.e. the condition A^ch 3> 1 is necessary. 

The method of Levitov and Shytov 2:] is based on clas- 
sic equations of motion for the electron density p{x, t) 
and current j(x,T) in imaginary time r = —it. In a case 
of wire one can write 

^ + t^^J{x,T), (11) 
OT OX 

J+D^-^-^ fd^'Pi^'^*)^^' (12) 
ox ox J £\X — X \ 

where a = e'^^/2Trh is effective one-dimensional conduc- 
tivity, e is the effective dielectric constant, and _D is a 
diffusion constant. The instanton is chosen in a form of 
a symmetric bounce, so that the source J7 in the continu- 
ity equation (|lip corresponds to the injection of one extra 
electron into the system at time ti = —itq at point x = 
—L/2 with its subsequent elimination at the same point 
at moment ^2 = itq: J\ = [(5(t-(-to) — 5(r — ro)](5(x-f-L/2). 

The crucial point of our approach is that we describe 
contunnelling through a diffusive wire by the same semi- 
classical equations ()11I12|) . but with modified source 

J2 = [<5(r + To)-5(T-ro)][,5(x + L/2)-<5(x-i/2)] (13) 

which describes simultaneous tunnelling of an electron 
and a hole via both contacts. 

The density p{x,t) should be real and even with re- 
spect to T — r, while the current j{x,T) should be 
purely imaginary and odd. p(x,t) and j{x,T) are de- 
fined on the interval — 1/2T < r < l/2r and obey peri- 
odic boundary conditions. Expanding p{x, r) and j{x, t) 
in Fourier series, we get p{x, t) — J2u ^) cos(a;T) and 
j{x,T) = j{x,LL!)sm(LUT), where the Matsubara fre- 
quency summation, as usual for Bose excitations, runs 
over even frequencies lu — 27rTn, with n = 0, 1, 2, . . .. 

If the wire is very long (Lq 00), the system of equa- 
tions (I11I12|) can be solved by the spatial Fourier trans- 
formation. Then, proceeding in the full analogy with [3], 
we obtain w 5^ exp{— S'2(T, V)}, where 



S2iT,V) = S2{T,T*)-2eVTS, 



dS2iT,To)/dTo 



= 2eV. 



(14) 
(15) 
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S2iT,To) 



E 



{uj + Dq^){uj + aq^Uq) 



(16) 



J2[uj, q) = -4i sin(wTo) sin(qL/2), (17) 



1 



U{q) = reg / —-dx = - In — . (18) 
J-oo e qa 

The semiclassical method, used above, is applicable, if 
^2 > 1- From (|16I17P it is clear, that 882/8x0 = for 
To = l/4r. Therefore we conclude that 



ro*(r,i/^o) = i/4r. 



(19) 



so that in the low-voltage case the summation runs only 
over the odd n = 2fc + 1: 



S2{T,V^{)) 



E 

fe=0 



dq 



16e'^Tsm^{qL/2)Ug 



27r (27rr(2/c + 1) + Dq^){2T:T{2k + I) + dq^Uq 



There are three different temperature ranges: T 3> ujn 

Wmin < r < Wmax; and T < tJniin, whcrC 



'^max 

-2 



aL-2[/(L-i) ^ (e/L)i?c(i), (20) 



ch- 



ln(L/a) 



<Wmax, (21) 



we consider these three ranges separately 

1. T 3> ti'max- Here the sum is dominated by fc ^ 1, 
u; - T; it can be shown that S2{T, ^ ^ 0) « 251 (T, V 
0), which means that in this temperature range the 
two-particle process looses a competition with the one- 
particle Coulomb anomaly one. 

2- Wmin ^ T <g; Wmax- Here the integral over q is 
dominated by g L^^, while the logarithmical sum over 
k is dominated by an interval < A; ^ Umax/T, so that 



S,{T,V^Q)^^\n{^) 



(22) 



and, with the help of ([^0]) . we arrive at the final expres- 
sion 

3. r ^ t^min- In this range presumably the elastic 
cotunneling should dominate. However, since Wmin 
Tloc, the localization effects may be important here, and 
we do not discuss this regime. 

In the case of high voltage one has Tq <^ 1/7", so that 
the infrared cutoff of the frequency summation in (|16p is 
I/to instead of T. As a result, we arrive at the same three 
regimes, but with substitution T ^ I/tq. In particular, 
in the most interesting regime 2: 

Q (rj. , 4L{eV) 

^2(T, To) ~ — 111 (^maxTo) , Tg = , (23) 

and we arrive at the high- voltage version of formula ^ . 

In the case of a finite wire the straightforward Fourier 
analysis of the problem is impossible, and one has to 



be more accurate. It is important that, for two-particle 
tunneling, the total charge of the wire is always zero, so 
that one does not have to take into account the effects of 
the total charge which otherwise would be important in 
a finite system at low temperature. This simplification 
allows one to follow the lines of the solution described 
above, with a substitution of the Fourier analysis by de- 
composition over a set of eigenfunctions of a certain linear 
operator. 

We present the action in a general form 

^-fEl-^+p^p: 



y 51 ^2 <! --g+ui/^vug + g+ug \ J2 = 



4sin^(ti;To) 



J2^+C/^2, (24) 



P = gj2, j ^avup^awugj2, g^ oj-dAu 



A = 



dx^ ' 



[Up]{x) = reg 



p{x')dx' 



Lo/2 



(25) 



The rules for regularization of the singular integral ([25]) 
are similar to those, given in 0. Within the logarithmic 
accuracy one can write [Up]{x) « U{x)p{x), where 



U{X) EE - 

e 



-min{(i/2-a;), A(a;)} 
{{L/2 + x), X{x)} 



2 . 
— mm 

a 



(26) 



and A(a;) ^ [d\n{p{x)) / dx]~^ is a characteristic scale of 
spatial variations of the function p{x). The formula pS)) 
is applicable, if the arguments of both logarithms are 
large. At A > io (EH) coincides with the result in [9j. 

Let us introduce a set of normalized right- 
eigenvectors Lpm{x) and eigenvalues Vm of the (non- 
Hermitean!) operator UA, obeying the boundary condi- 
tions d(p/dx\,-^^j.]^^^2 — Oj corresponding to a vanishing 
electric field at the ends of the wire. The equation and 
the normalization condition for the eigenfunctions read 



Lo/2 



Aifimix) = VmU Vm(a:^), 

ipm{x)U^^iPm'{x)dx = 5mm'- 



(27) 

(28) 



-La/2 



There is only one zero eigenvalue wq = 0, the correspond- 
ing eigenniode ipo{x) = const describes the equipotential 
distribution of charge in the wire. For the nonzero modes 
(m =1,2,...) in the leading logarithmic approximation 



^Prn [x) 



l2UmiO) 

Lo 

Vm = 



COS 



nm mnx 

2™2 



n'm'Um{0)/Ll 



(29) 
(30) 
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where Um{x) is given by ((26)) with A(a;) = Lo/{m + 1) 
Then, for Hermitean operator ■ U, one has 



LO — av,. 



(31) 



so that the action can now be rewritten in the form 



(32) 



/T NO /2C^m(0) . 7rm . nmL 

(J2 • ^m) = 2W— ^— sm — sm-^. (33) 

Thus, only odd to = 2p+ 1, p = 0, 1, . . . are relevant, and 



IGe^sin^wTo ;72p+i(0)sin 



2 7r(p+l/2)L 



LqU! 



p=0 



7r^(2p+l)^o-i72p+i(0) 



(34) 



Since ?7m(0) depends on to logarithmically, one can re- 
place Um{0) U = UmiO). Here to is the characteristic 
value of TO (or of 2p + 1), corresponding to those terms 
in (j34p . that give the principal contribution to the sum. 
The value of to will be found a posteriori. As a result 



^2 



16e^t/ sin^ wtq 



in 



2 7r(p+l/2)L 



7r^(2p+l)^art/ 

p=0 ^ 



(35) 



For F ^ we have again Tq = 1/4T. In the range of our 
interest (wmin ^ T ^ ij^max) one can neglect the term tu 
in the last denominator of (1351) and write 



wire has no effect on the process. This is not the case 
for one-particle tunneling, where the charge tends to pro- 
liferate equally in both directions, and the size-effect is 
important. These problems will be discussed elsewhere. 

For interpretation of the result ([6]) let us think of the 
relevant stretch of the wire between the two contacts as 
a sequence oi N — L/S^ grains - small pieces of length ^, 
each, connected by conductances gcs ~ 1- Strictly speak- 
ing, the perturbation theory is only applicable, if g^s ^ 
1; however, for g^s 1 it should still give qualitatively 
correct estimates. Then, using the results of [lOj, we esti- 
mate the effective conductance of this string of grains as 

Geff g'{g.sN'TyEciOY - g\LT/S,Ec{S,)f^'^. 
which roughly agrees with ([5]). 

In conclusion, we have shown that low-temperature dif- 
fusive transport through a wire, connected to the leads by 

two tunnel junctions, is realized as an inelastic cotunnel- 

(2) 

ing process, the effective conductance G\g obeys a power 
law with the L-dependent index. Our result ^ seems to 
be in agreement with the data[lj. In contrast with the 
baUistic case, where a <C 1 for A^ch ^ 1, our a* can 
be quite large. The expression ^ for a* is apparently 
similar to the formula a = Z(0)/(/i/2e^), used in the 
phenomenological " environmental theory" [O, [i3| ■ The 
latter approach, however, was not able to reproduce the 
correct exponential behavior ([3]) for the single-junction 
setup; neither can it produce specific value R{L) for the 
effective impedance Z{Q). 
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(36) 



Analyzing the series in ([M]) . we find to = 2p + 1 ~ Lq/L. 
Hence 



U ' 



e a 



2^2 In k 

a_ 



{i/L)Ec{L). (37) 



The high- voltage modification of (|36p is obtained exactly 
in the same way, as it was done for the case of infinite 
wire. 

Thus, in the leading logarithmic approximation, the 
finite-size effects do not modify the result (|6]). The rea- 
son is in the special geometry, characteristic for the two- 
particle cotunneling process: The initial dipole-like dis- 
tribution of charge shrinks, so that the charges always 
move towards the center of the wire, and not in the op- 
posite direction; therefore the presence of the ends of the 
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